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Abstract. In this article we study the large-time behavior of the solution to a general Rosenau 
type approximation to the heat equation [15], by showing that the solution to this approximation 
approaches the fundamental solution of the heat equation at a sub-optimal rate. The result is valid 
in particular for the central differences scheme approximation of the heat equation, a property which 
to our knowledge has never been observed before. 
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1. Introduction 



r ^~l ■ 1.1. Rosenau Regularization of the Chapman-Enskog Expansion. In his seminal work [IB] . 
!> | Rosenau has proposed a regularized version of the Chapman-Enskog expansion of hydrodynamics. 
^ • This regularized expansion resembles the usual Navier-Stokes viscosity terms at low wave-numbers, 
. but unlike the latter, it has the advantage of being a bounded macroscopic approximation to the 
linearized collision operator. The model is given by the scalar equation 

i v 



df gg(/) 
v ' ' at dv 



f(0 



1 + e 2 m 2 ^ 2 



where /(£) denotes the Fourier transform of f(v), while /(£) V denotes the inverse Fourier transform. 

The operator on the right hand side looks like the usual viscosity term ef vv at low wave-numbers 
£, while for higher wave numbers it is intended to model a bounded approximation of a linearized 
collision operator, thereby avoiding the artificial instabilities that occur when the Chapman-Enskog 
expansion for such an operator is truncated after a finite number of terms. 
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One of the advantages of this regularization is that it is equivalent in small frequencies to the heat 
equation (and then regularizes initial data), whereas high frequencies behave as an absorption term. 
In particular, information will travel at finite speed, and tools from hyperbolic equations can be used. 
Note that the right side of can be written as 



(1.2) 



e£ 2 

s /(£) 



-\ v 



1 + e 2 m 2 £ 



(em) 



1 m - m 



1 + e 2 m 2 ^ 2 



— [M e -*f-f] 
me 



where e = me, * denotes convolution and 

(1.3) M y («) = ^e-M/T 

is a non-negative function satisfying ||M 7 ||^i = 1. In other words, the Rosenau approximation 
consists in substituting the linear diffusion equation 

d.4) ! ( „ ()= ^ M) 

with the linear kinetic equation 

(1-5) ^(v,t) = ^[M e *g(v,t)-g(v,t)] 

in which the "Maxwellian" M e [7] is given by (|1.3|) , 

Equation (|1.1|) has then been studied by S. Schochet and E. Tadmor in [T7] in the context of 
vanishing viscosity solutions of scalar conservation laws. H. Liu and E. Tadmor then studied in 
|13] the Burgers equation with a Rosenau-like nonlocal viscosity term and proved the existence of 
a so-called critical threshold, namely the existence of a critical value of the total variation of the 
initial condition triggering or not a finite-time blow-up of the solutions. This nonlocal term was also 
used by C. Rohde to model capillarity effects close to fluid- vapour phase transitions in Navier-Stokes 
equations [15] . Interestingly, this non-local approach seems more physically relevant than the original 
Korteweg's idea to use the heat operator for describing such phenomenons. 

More recently, it has been replaced in a general context of system of hyperbolic balance laws with 
non-local source term by R. Colombo and G. Guerra in the two companion papers OS]. They showed 
in particular the well posedness in L 1 (globally in time) of the Cauchy problem for (jl.l|) together 
with uniform stability estimates, in the framework of viscosity solutions. Finally, C. Besse and T. 
Goudon obtained recently in [3] a system of macroscopic equations with a nonlocal Rosenau term by 
computing the diffusion limit of a linear space inhomogeneous kinetic equation. A similar approach 
was also used by T. Goudon and M. Parisot in [T2] on a kinetic equation describing the behavior of 
a multi-species gas of charged particles to obtain a macroscopic model with a Rosenau-like nonlocal 
term. 

Despite the previous studies, it is not completely clear and well established from a mathematical 
point of view if the correction proposed by Rosenau is a good approximation to the linear heat 
equation. In particular, it is not clear whether or not the large-time behavior of the solution to the 
Rosenau approximation agrees with the large-time behavior of the linear diffusion equation. 

The aim of this article is to give an answer to the previous question, and to underline that the 
Rosenau approximation can be viewed as a particular case of a general approximation to the heat 
equation by means of a linear kinetic equation of type (|1.5|) . provided the background density M £ is 
a probability density function of zero mean and variance proportional to e 2 . In particular, it will 
be shown that, in a certain metric equivalent to the weak*-convergence of measures, the distance 
between the solution to the heat equation and the solution to the kinetic equation can be bounded 
uniformly in terms of e and t, provided the background density has a sufficiently high number of 
moments (typically more than two). 

The plan on the article is as follows. In Section [21 using tools of the kinetic theory of rarefied gases, 
we will introduce a possible derivation and the main features of the Rosenau approximation with a 
general kernel. This allows to describe well known models such as the central differences schemes 
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for the heat equation. In Section^ this kinetic formulation will be used to obtain explicit solutions 
to the Rosenau equation (|1.5|) . using both Fourier transform and Wild sums. Last, in Section [U 
we will investigate the large time behavior of the solutions to (|1.5|) . showing that the convergence 
towards the fundamental solution occurs, in a suitable Fourier-based metric (see next subsection), at 
a suboptimal rate (compared to the heat equation). Finally, we will combine these results to show 
that strong convergence in L 1 towards the fundamental solution to the heat equation is obtained 
after a suitable regularization of the Rosenau equation, obtained by discarding its singular part. 

1.2. Functional Framework. Before entering into the main topic of this paper, we list below the 
various functional spaces used in the following. For p £ [1, +oo) and q E [1, +oo), we denote by IP q 
the weighted Lebesgue spaces 

L p q := j/ : R R measurable; ||/||^ := J \ f(v)\ p (1 + v 2 ) q/2 dv < ooj . 

In particular, the usual Lebesgue spaces are given by 

IF := Ll 

Moreover, for / £ ij, we can define for any k < q the k th order moment of / as the quantity 

m k {f) ■= / f(v) \v\ k dv < oo. 



For s S N, we denote by W s ' p the Sobolev spaces 



P dv < oo 



W*:=\feL'; \\f\\ p WSjp := £ / |/ (fc) 

I |fe|<s 

If p = 2 we set H s := W s ' 2 . 

Given a probability density /, we define its Fourier transform J- V (f) by 

Mf)(0 = AO ■= I e~^ v f(v) dv, V£ G R. 
The Sobolev space H s can equivalently be defined for any s > by the norm 

:=II^(/)IL| • 

Is 

The homogeneous Sobolev space H s is then defined by the homogeneous norm 

11/6. - / \t\ 2s m 2 

Finally, we introduce a family of Fourier based metrics in the following way: given s > and two 
probability distributions /i and /2, their Fourier based distance d s (/i>/2) is the quantity 

7i(6-M£) 



d s {h,f2j ■= sup 

5eR\o l?r 

This distance is finite, provided that /i and /2 have the same moments up to order [s], where, if 
s ^ N, [s] denotes the entire part of s, or up to order s — 1 if s £ N. Moreover d s is an ideal metric 
[B]. Its main properties are the following 

(1) For all probability distributions /i, f^, f-s, 

ds(fi*h,f2*f3)<d s (f l ,f 2 ); 

(2) Define for a given nonnegative constant a the dilatation 



Then for all probability distributions fx, /2, and any nonnegative constant a 

d s (h, a , h,a) = a s d s (fi, / 2 ). 
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The c^-metric is related to other more known metrics of large use in probability theory [TU] . 

2. A Kinetic Description of the Rosenau Approximation 

Let v € M denote velocity, and let us assume to have at time r > a system of particles immersed 
in a background. Let us suppose that the number of particles in the system is sufficiently large to be 
studied by means of statistical mechanics, namely by giving the velocity distribution f(v,r) at time 
t > 0. Moreover, let us assume that the main phenomenon which can modify particle's velocity is 
the interaction of particles with the background. Let M(w) denote the fixed in time (probability) 
distribution of the particles of the background, which we will assume of finite variance 7 2 . Assume 
that the collision process of a particle with velocity v with a background particle with a velocity w 
generates a post-collision velocity v* given by 

(2.1) v * = v + w . 

Then, in a suitable scaling [7], the effect of interactions (|2,ip on the time- variation of the density 
f(v,r) can be quantitatively described by a linear Boltzmann-type equation, in which the variation 
of the density is due to a balance between a gain and loss terms. 

More precisely, for a given number v, we take into account all the interactions of type (|2.1|) which 
end up with the number v* (gain term) as well as all the interactions which, starting from the number 
v, lose this value after interaction (loss term). The balance equation for the density of particles can 
be fruitfully written in weak form. It corresponds to say that the aforementioned interaction process 
on particles modifies the solution f(v,r) according to 

(2.2) — / ip(v)f(v,T)dv = \ [ (<p(v*)-<p(v))f(v,T)M(w)dvdw, 
dr Jm jr 2 

where the constant A > denotes the intensity of the variation process, and <p(v) is a smooth 
function. Note that choosing (p(v) = 1 shows that, independently of the background distribution, 
f(v,r) remains a probability density if it so initially 

f(v,r)dv = / fo(v)dv = 1. 
Jm 

This is in general the unique conservation law associated to equation (|2.2|) . 

From now on, let us assume in addition that the probability distribution of the background is 
centered, and its variance depends on a small parameter e > 0. To emphasize this dependence, we 
will denote this distribution by M £ (w) = e M(e~ l vj). Then M £ (w) satisfies 

(2.3) / M £ (w) dw = 1, f w M £ (w) dw = 0, / w 2 M £ (w)dw = e 2 . 



The weak formulation (|2.2|) yields immediately the time evolution of the moments of /. Taking 
ip(v) = v one obtains 



d_ 

~d7 

Moreover, if <p(v) = v 2 



vf(v,T)dv = X / w M £ (w) dw / f(v,r)dv = 0. 



d 



(2.4) — / v z f(v,T)dv = X / (v*y-v^ f(v,T)M £ (w)dvdw = \e 1 ^ 1 . 

dr Jm. Jm? 1 J 

Thus, the second moment of / grows linearly with respect to time and depends on e. One way to 
avoid this dependency is to scale the time properly. Setting t = e 2 t and introducing a new particles 
distribution function g £ such that g £ (v,t) = f(v,r) gives according to (|2.4|) . 



jt L y2aM t)dv = hiL v2f(v > r ) dv = A7 2 , 

and the second moment of g £ does not depend on e. 
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The distribution g £ is a weak solution to 

(2.5) / g £ (v,t)<p(v)dv = \ [ [<p(v + w) - tp(v)] g £ (v,t)M £ (w) dv dw, 
dt Jm. £ z Jm? 

with g £ (v,t = 0) = go(v) = fo(v). Note that, since 

ip(v + w)g £ (v, t) M e (w) dv dw = / tp(z) g £ (v,t) M £ (z — v) dv dz 

2 

<p(z)M £ *g £ {z) dz, 

JR 

equation ()2.5|) can be rewritten as 

d f A f 

(2.6) — / g £ (v,t)tp(v)dv = -s / <p(y) (M £ * g £ (v) - g £ (v)) dv, 
dt jr e z jr 

which is the weak form of (|1.2|) . Hence, if the background distribution is given by (|1.3|) . equation 
(|2.5|) coincides with Rosenau's approximation (|1.5p . where cr 2 = A. Of course, other choices of the 
background are possible, and, provided conditions (|2.3|) are satisfied, the evolution of the moments 
of the solution to (|2.5|) up to the second order do not depend on the background distribution. 

It is then easy to see that, in the case in which M £ coincides with the "Maxwellian" (|1.3p . g £ is an 
approximate (at the leading order in e) weak solution of the heat equation (|1.4|) . Indeed, given that 
e is small enough, one can Taylor expand f(v*), where v* = v + w, to obtain 

tp(v*) = ip(v) + w <p'(v) + ^- tp"(v) + ^(pW(v)w 3 
for v G (u, u + w). Using this relation in (j2.5|) . one obtains 

(2.7) i- / ^(i;,*)^)^ = ^ / 5 e («,i)^ / («)^ + J R(e) 

ai jr 2 

where the remainder satisfies 

(2.8) \R(e)\ < A||^(3)|| Loo J_ f w *M e (w)dw = eh\tp®\\ Jja o [ w 3 M(w)dw. 

3! e z Jr 3! Jr 

Remark 1 . One can easily notice thanks to (|2.7p that equation (|2.5|) is an approximation to the 
heat equation (|1.4|) (with diffusion coefficient a 2 = Xj 2 /2), provided that the remainder converges to 
zero as e — > 0. In order to have this convergence, it is enough that the distribution of the background 
is such that some moment of order greater than two remains bounded. We will use this in the 
following to obtain various approximations of the heat equation, just changing the distribution of the 
background. 



Remark 2. The Rosenau type kinetic equation (|2.5p is such that mass, momentum and energy of its 
solution have the same evolution of the corresponding moments of the solution to the heat equation 
(|1.4|) . A further interesting analogy with the heat equation is given by studying the evolution of 
convex functionals along the solution. Let 3>(r), r > be a (regular) convex function of r. Then, 
using equation (|2.6|) we obtain 

d f Hg £ (v,t))dv= [ ^'(g £ (v,t))^4^dv 



dt Jr Jr dt 



4 / &(g E (v,t))(M £ *g £ (v)-g £ (v)) dv. 



e 2 



Thanks to the convexity of $(•), for r, s > 

$'(s)(r -s) < $(r) - 

and one obtains 

/ ®(g £ (v,t))dv<± [ ($(M e * g £ (v)) — <&(g £ (v))) dv. 
dt Jr e A Jr 
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Now, use the fact that M £ is a probability distribution, so that by Jensen's inequality 

$(M £ *g e (v))dv = /$( / g £ (v - w))M £ (w)dw j dv 
jr \Ju J 



< / <Z>(g £ (v - w))M e (w) dv dw = / $(g e (u))du, 

JR 2 JR 

and one concludes with 

d f 

- / $(g e (y,t))dv<0. 
at Jr 

Thus any convex functional is non-increasing along the solution to the Rosenau type kinetic equation 
(jl.5p . in agreement with the analogous property of the heat equation. 



Remark 3 . A leading example of Rosenau type approximation is obtained by assuming that the 
background distribution M £ is a balanced Bernoulli distribution 

(2.9) M £ {v) = 1[S (v + e 1 ) + S (v-e 1 )}, 

where as usual 5q denotes a Dirac mass concentrated at v = 0. Note that the "Maxwellian" (|2,9p 
satisfies (|2.3p . Then, according to (|2.6p . the distribution g £ solves 

(2.10) -^(v,t) = ^[g £ (v - e 7 ,t) - 2g £ (v,t) + g £ (v + e 7 ,t)] . 

Now, let us fix A = 2 and 7 = 0". Then, given a parameter Av > and a uniform grid = 
iAv, i £ Z on R, by setting #j(t) := g £ (vi,t) and e = Av/a, equation ()2.10p reads 

2 5i+i(*) - 2 9i(t) + 9i-i(t) 



CI' 



Av 2 

This is exactly the classical semi-explicit second order central differences scheme for the heat equation 
(|1.4p . It follows that this scheme furnishes an approximation to the heat equation which satisfies 
all the properties outlined in Remarks [T] and \2\ as well as the weak convergence towards the exact 
solution, which follows easily by applying the estimate 



Remark 4 ■ The previous example introduces into the matter a background with a non regular 
distribution (presence of point masses). This situation clearly differs from the standard Rosenau 
approximation, characterized by the regular density (|1.3|) . Due to its wide applications in numerical 
simulations, the study of the properties of the Rosenau approximation given by this type of back- 
ground turns out to be important. For this reason, in addition to the study of the classical Rosenau 
model, in what follows we will deal also with the Bernoulli type background, outlining when possible 
the principal differences. 

3. Representations of the Solutions to the Rosenau Equation 

The Rosenau approximation (|1.5p of the heat equation is a linear kinetic equation of Boltzmann 
type. Therefore, we can resort both to linear kinetic theory and to the theory of linear diffusion 
equations to obtain explicit representations of the solution. We shall present here two equivalent 
ways to construct these solutions, each one giving a different insight on the behavior of g £ . In the 
rest of this Section, to avoid inessential difficulties, we will assume that the initial value go{v) is a 
probability density with finite moments up to a given order (in general more than two). 

3.1. Wild Sums. The first method we introduce to obtain an explicit representation of the solution 
g £ is closely related to the so-called theory of Wild sum expansions of the Boltzmann equation. Wild 
indeed proved in |19j that one may represent solutions to the nonlinear Boltzmann equation for 
Maxwell molecules by using convergent power series. This method is particularly simple when one 
deals with linear convolution equations such as (|1.5|) . To start with, let us set 

h £ (v,t) := exp g £ {v,t). 
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Then the Cauchy problem for (|1,5|) can be rewritten as a fixed point problem (which is nothing but 
a Duhamel formula). Let h —> <& £ (h) define the map 

A C l 

(3.1) <5> e {h) = g + — M £ * h{s) ds. 

e A Jo 

Then g £ (y,t) solves (|1.5p exactly when $> e (h e ) = h £ . 

To find the fixed point, it is then sufficient to make a Picard iteration. Starting from h £ := go 
one defines, for n E N 

(3.2) h e n+1) = ®e (h e n) ) . 
By recurrence, for n > 1 it yields 

(3.3) hW(v,t) = h^Hv^ + i^j -M^*g (v), 



where we use the shorthand M* n := M * ■ ■ ■ * M (n times). Clearly, for n > 

hi n+1) (v,t)-hi n \v,t) >0. 
Hence the sequence (h e n \v, t)] , being bounded in L 1 , converges towards h £ (v,t) > when 

V / n>0 

n — > oo, where according to ()3.3p 



/ A t \ n 1 

h £ (v, t) = go(v) + E — "I M * n * 9o(v). 



n>l 



nl 



By passing to the limit n — > oo in (|3.2|) . one obtains that h £ is a (nonnegative) fixed point for 
This procedure allows to write the (nonnegative) solution g £ to the Cauchy problem for (|1.5|) as 

(3.4) g £ (v,t) = e- xt ^g (v) + e- xt ^ ^ f^Yl-M^*g (v). 

n>l 



3.2. Fourier Transform. The standard approach to solve the heat equation on the whole space is 
to use Fourier transform. This is even more obvious for equation (|1.5p . Indeed, passing to Fourier 
variables in this equation shows that the Fourier transform of g £ is given by 

(3.5) &(t,t)=go(£)exp(-MOt), V(e,t)EMxR+, 

where is the Fourier transform of the initial density and 

i - K(0 



(3.6) A E (0 := A 



£ 2 



Since M £ is a probability density, |M e (£)| < 1 for all £ E K, which implies that the real part of A £ is 
nonnegative. Moreover, thanks to condition (|2.3p . if M £ possesses more than two moments bounded, 
one can write 

where the o(£ 2 ) term is considered for £ — > (which will be the case until the end of the article). 
Finally, it is possible to invert the Fourier transform (|3.5|) and obtain in a distributional sense 

(3.7) g £ (v,t)=g *P £ (;t)(v), 

where we set 

P £ (v,t) =exp{-A e (Ot} V (v). 
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Figure 3.1. Schematic representation of the convolution steps of the central differ- 
ence kernel. Each node is a Dirac mass, weighted by binomial coefficients. 

Clearly, representations (|3.4|) and (|3.7|) have to coincide. Indeed, applying a Fourier transform to 
Q32D yields 

■\t\ n 1 



(3.8) 



-Xt/e' 



-Xt/e 2 



5o(0 + e 



-Xt/e 2 



E 

n>l 



n! 



3 M £ (C)At/£ 2 



1+1 



exp(-A £ (£)i) go (0. 



where ^4 £ is given by (|3.6p . 



Remark 5. According to relation (|3.8p . any Rosenau-type approximation to the heat equation is 
such that the smoothness of the solution is lost, independently of the regularity of the kernel M £ . 
Indeed, the fundamental solution to equation (|1.5p (namely the solution to the Cauchy problem 
obtained from the initial value go = 5q) can be represented in Fourier variable by 



Xt/e 2 M E (i)Xt/e 2 



1 +e 



-Xt/e 2 



(3.9) =:G 1 (£,t) + G 2 (t). 

If M £ £ L 1 , then according to Riemann-Lebesgue Lemma, M £ (£) — > when |^| — > oo and then 
|Gi(£,i)| — > 0. The loss of smoothness comes from G-2,{t) which is the Fourier transform of a Dirac 
mass (of weight depending on time). According to (|3.9p . the singularity disappears exponentially 
when e — > or time goes to infinity. 

If one considers the Rosenau kernel (|1.3|) with A = a 2 , 
and for all t, e > 0, 

Thus Gi(,£,t) decreases exponentially fast towards zero when |£| 
comes only from Gi. 

On the contrary, if one considers the central differences kernel 
M £ (£,) = cos(ea^) and for all t, e > 0, 



oo, and the loss of smoothness 
with A = 2 and a > 0, then 



Gi(£,t) =e 



-2t/e 2 



exp [ — cos(eo"^J 
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which does not converge at infinity. Then, if the kernel has some singularities, the solution to the 
Rosenau equation has at most the regularity of the initial condition. This enlightens a substantial 
difference of the Rosenau approximation with respect to the heat equation, in which the diffusion 
operator is such that the solution is instantaneously smoothed out. The representation (|3.4|) can be 
used to understand the behavior in time of the solution obtained from the central differences kernel. 
Indeed, each convolution by M e consists in creating a weighted sum of each point (i.e. Dirac masses) 
of the grid (see Figure 13. lj) : 

M £ * Sq(v) = - (S (v - ea) + S (v + ea)) , 

M* 2 * 5 (v) = ~ 5 Q (v - 2ea) + ~ 5 (v) + J S (v -2ea), 

1 3 3 1 

M* 3 *8 (v) = -5 (v - 3ea) + - 5 (v - ea) + - 5 (v + ea) + - 5 (v + 3ea) 

o o o o 

and so on. Thus, the fundamental solution of the central difference scheme is given by 

9e(v,t) = e-^ 2 5 (v)+e-^ 2 j:( 2 4) n ^9i n \v), 

n>l V£ J U 



where, for n G N we defined 

9? n) (v) = Q)"e(J) Mv + 2kea)+5 (v-2kea)}, 

gf n+1 \v) = (\y +1 J2 Qlll) [5 (v + (2k + l)ea) + 5 (v - (2k + l)ea)] . 

4. Asymptotic Behavior of Solutions to the Rosenau Equation 

In Section [5] the similarities between the Rosenau-type kinetic equation and the linear diffusion 
equation have been enlightened. In particular, as well known from the classical literature on numerical 
approximation of partial differential equations, equation (|1.5|) provides a consistent approximation 
of the heat equation in a fixed time interval. It is not known, however, if equation (|1.5p still realizes 
a good approximation for large times, and, in case of a positive answer, in which way the difference 
between the solutions of the Rosenau and diffusion equations can be estimated with respect both to 
time and the parameter e. 

In this Section, we will furnish a partial answer to this question. To clarify our intent, we will 
briefly resume various well-known facts about the large-time behavior of the solution to the heat 
equation posed in the whole space. Let g = g(v, t) be a solution to the heat equation (|1.4|) and define 
the heat kernel £l a by 

(4.1) n„(v, t) := ^L= exp (--jQ , V(«, t) E R x R + . 

This fundamental solution of (|1.4|) represents an intermediate asymptotics of a large class of solutions 
to the heat equation. The recent review article [2] gives a precise state of the art on this topic. To 
make this concept more precise, we define the Boltzmann entropy of / as 

«(/):= / f(v) log f(v)dv. 



Then it can be shown (see e.g. [TB]) that g(v,t) behaves as the heat kernel when t — > oo, provided 
that the initial condition go is of finite kinetic energy and entropy: 

go E L\ and W(g ) < oo. 

Moreover, the rate of convergence towards the fundamental solution can be computed in L 1 norm 

(4.2) ii^-n^)^ <_^, 
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where C is an explicit constant. The bound (|4.2|) is sharp. A marked improvement of the constant 
in (|4,2|) has been recently obtained in pQ, by selecting well parametrized Gaussian functions, char- 
acterized either by mass centering or by fixing the second moments or the covariance matrix of the 
solution. 

Most of these results follow by considering the equivalence between the heat equation and the 
linear Fokker-Planck equation 

which has a stationary solution given by the Gaussian 

(4.4) Lo a (v) :=n a {v,t = l). 
Setting T(t) := e 2 '" 1 / 2 , 

we can in fact consider the change of variables 

(4.5) u(v,t) :=Jg(Jv,T(t)). 

Then u(v, t) is solution to the Fokker-Planck equation (|4.3|) as soon as g is a solution to the heat 
equation (|1.4|) . The converse also holds true. It has been shown in [5] that the solution to the linear 
Fokker-Planck equation converges for large times to the stationary solution u) a and the following 
asymptotic behavior holds 

\\u(t) -uvlUi < Ce'\ Vt > 0. 

Using the self-similar scaling (|4.5|) . it is clear that this result is equivalent to the rate (|4.2p . and 
once more, the bound is sharp. We shall use a similar rescaling in our study of the intermediate 
asymptotics of equation (|1.5|) , 

Our goal in the following is to prove that the same kind of result holds for the Rosenau-type 
approximation to the heat equation. Due to the generality of this approximation, which depends 
of the background distribution M, we will obtain a weaker convergence (with respect to the L 1 
distance) to the fundamental solution, as discussed in Section [21 Remark [5l The L 1 distance will be 
here substituted by the suitable weaker Fourier based metrics d s described in Section 11.21 which are 
particularly adapted to the convolution structure of the kinetic equation (|2.5|) . This family of metrics 
has been introduced in the article [TU] to study the trend to equilibrium of solutions to the space 
homogeneous Boltzmann equation for Maxwell molecules, and subsequently applied to a variety of 
problems related to kinetic models of Maxwell type. For a more detailed description, we address the 
interested reader to the recent lecture notes [BJ. 

These metrics can be easily and fruitfully applied to the study of the large-time behavior of the 
heat equation. Indeed, we have seen that in Fourier variables, the solution to the heat equation (|1.4|) 
with initial condition go is 

9&t)=9b(0 exp(-cT 2 £ 2 t) . 

Thus, if go is of finite mass, the solution will converge pointwise towards when t — > oo. A way to 
consider a nontrivial limit distribution is to make the change of variable £ —¥ V(t) £ with V(t) := 
(1 + 1) -1 / 2 . Then, the scaled distribution will converge towards the Gaussian distribution 



g (0) exp (-v 2 t 2 ) , 

and one can improve the rate (|4.2|) in Fourier distance, as soon as the initial datum possesses finite 
moments of certain order equal to those of the Gaussian [11] . From now on, let us set 

(4.6) u5 ff (£) = exp {-cr 2 e) , 

the Fourier transform of the Gaussian distribution (|4.4|) . We have the following result. 

Proposition 4.1. Let V(t) := (1 + 1) -1 / 2 . For a given s > 1, let < go G L\ a distribution of unit 
mass. Then, if g(v,t) is the unique solution to the Cauchy problem (|1.4p with initial condition go, 
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the scaled distribution h(v, t) := V(t) 1 g (V(t) 1 v, t) is such that, for all t > 

1 



(4.7) 



d s (h{t),u a ) < 



; d s (g ,uj a ). 



(i + ty/ 2 

Proof. Thanks to the scaling properties of the Fourier transform one has 

\g(V(t)Z,t)-uT a (0\ 



d s (h(t),u} a ) = sup 
£GK\0 

= sup 



go 



(y(t)O-exp {-a 2 eV 2 {t) 



< V(t) s sup \m) ~ exp (-° 2 f) 



1 



exp I — a 2 ^ 2 



l+t 



□ 



Remark 6 . Taking a Dirac mass as the initial condition in the Cauchy problem for equation 11.41 
shows that these rates are optimal. 

4.1. Approximate Solutions with Finite Energy. To recover the asymptotic behavior of the 
solution to the Rosenau-type approximation (|2.6|K we are going to apply a technique similar to that 
used in Proposition 14.11 Let us consider again the scaling £ -> V(t)£ with V(t) := (1 + t)~ 1/2 , and 
define 

h(v, t) := Vit)- 1 g (Vit)- 1 v, t) and h £ (v, t) := V(ty l g £ (y(t) -1 v, t) . 
We choose in this section an initial datum < go £ L\ such that 

go(v)dv = l, / vgo(v)dv = 0, / v 2 go(v) dv = E < oo , 



(4.8) 

and A7 2 /2 = a 2 . Then the quantity c?2 (goikv) is bounded, and, according to (|4.7p one obtains 



(4.9) 



d 2 (M*W) < ^7 d 2 (5o,^) + d 2 (h E (t),h(t)),Vt>0. 



Hence, we have to estimate the Fourier distance between the exact scaled solution h and its approxi- 
mate counterpart h e . We shall use for this the representation (|3.5p of g e in Fourier variables. Let us 
introduce a positive parameter R. Since the mass of the initial datum go is equal to 1, one has 



1 



d 2 {h £ (t),h(t)) < sup 
2 



exp <^ -A 



< ^ + SUp 

2 1 

< -^j + sup 72 



exp <! 



exp I —a 2 ^ 2 



l + t 



i? 2 



I£I<R 



c 2 



exp ^ — cr 2 £ 2 



f }> - exp i -a 2 f 



t 



l + t 
a 2 e 



A, 



VT+t 



Using the elementary inequality |1 — e x \ < \x\, one finally has 



l + t 
- 1 



(4.10) 
where 



d 2 (h £ (t),h(t)) < — + sup D e (£,t), 



A 1 - M F 



\Z\<R 
VT+t 



l + t 



In order to compute this last quantity and to obtain the rate of convergence, let us specify both 
the background distribution and the constants 7 and A. 
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Central Differences. As in Remark [3\ let us set A = 2 and 

M £ {v) = -[5o(y + eg) + S (v - ea)]. 
Then M e (£) = cos(eo"£), and thanks to Taylor's theorem, there exists 9 6 (0, 1) such that 



(4.11) 



t eal; 



3(l + t) 3 / 2 



5111 



1 

y/T+t 



< 



Gathering inequalities (|4.1U|) and (|4.11|) yields for R > the inequality 



d 2 {h £ {t),h{t)) < 



+ 



(eoRf 



R 2 (1 + t) 2 3 



It just remains to optimize on R to obtain 
(4.12) 



d 2 (h £ (t),h(t))<J^e v ' 



1 + t 



Inequality ()4.12p shows that the second term in (|4.9p has a slower rate of decay with respect to the 
first one. In contrast to the exact solution h, which is known to converge in d 2 towards the self-similar 
one at rate (1 + 1) -1 [11], we found that the approximate solution converges towards the exact one 
at rate (1 + t)~ 1 / 2 . We just stated the following 



Proposition 4.2. Let < go £ L\ satisfy (|4.8|) . If g £ is solution to the associated Cauchy problem 
for the Rosenau-type approximation \2. b)) with the central difference kernel 112. 9\) . then the scaled 
distribution 

h £ (t,v) := Vit)" 1 g £ (y(t)~ l v,t)j 
for V(t) := (1 + t)~ 1 / 2 verifies for all t > 



d 2 (h e (t),u) a ) < 



3a 2 V~t 
- £ Z 2 ( 50 , W(7 ) + W_ e __. 



1 + t 

Rosenau Regularization. In this case the Fourier transformed "Maxwellian" reads 

1 



M £ (0 



and A = a 2 . We obtain 



(ea) 2 t 
1 + t 



l + (ea) 2 e 7 

e 



l + t + (ea) 2 e 



(eo-Rf 



Inserting this expression in inequality (|4.1U|) gives for R > 

It just remains to optimize on R to obtain the rate of convergence in d 2 norm of h £ towards h: 

+ Vi 



(4.13) 



la" 



d 2 (h £ (t),h(t)) < \lY £ \+ r 



Hence, this approximation has exactly the same order of convergence than the one given by the 
central differences kernel. We proved 

Proposition 4.3. Let < go G L\ satisfy (|4.8|) . If g £ is solution to the associated Cauchy problem 
for the Rosenau-type approximation \2. b}) with the Rosenau kernel then the scaled distribution 



h £ (t,v) := V{t)- 1 g s {y{t)- 1 v l t) 
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for V(t) := (1 + 1)- 1 / 2 verifies, for all t > 



1 a 2 \ft 

d 2 {h £ (t),u) a ) < -——d 2 {go,uJa) + \Hr^- 



1+t v ' V 2 1+t 

Remark 7. Note that as a by-product of our analysis of the large-time behavior of the solution to 
the Rosenau-type approximation of the diffusion equation we proved that the c?2-distance between the 
solution of the diffusion equation and its numerical approximation by the central difference scheme 
is uniformly bounded in time with respect to the parameter e, as given by (|4.12j) . 

4.2. Approximate Solution with Finite Fourth Order Moment. In Subsection 14.11 the con- 
vergence results required the boundedness of the ^-distance between the initial value of the Rosenau 
equation and the Gaussian density. This boundedness has been achieved by assuming that the initial 
datum is of finite energy. Indeed, moments play an important role with respect to convergence. We 
will improve here the results of Subsection l4.ll provided that the initial condition has finite moments 
up to order four. We will see that this allows to compute the rate of convergence (in time) towards 
self-similarity with very few assumptions on the kernel M £ . 

Without loss of generality (thanks to scaling and translational invariance), we can choose the 
initial distribution < go S L\ to satisfy 

(4.14) / g (v)cp(v)dv = (1,0,1,0,/x), 



where ip(v) = (1, v, v 2 /2, v 3 , v 4 ), and fj, is a positive constant. If g is solution to the Cauchy problem 
for equation (|1.4|) with initial condition go, then on the one hand 

(4.15) / g(v,t)<p(v)dv = (1,0,1 + to- 2 , 0, /i + 24t + Y2a 2 t 2 ). 

Jr 

On the other hand, a solution g e to equation (|2,6|) with same initial datum verifies 

(4.16) / g £ (v, t) (p(v) dv = (1,0, 1 + ta 2 , 0, /x + (24 + B £ )t + UaH 2 )), 
Jr 

provided that M e follows assumption (|2.3|) and where 

B e :=4 / M £ (v)v 4 dv. 
e z Jr 

We can see that the fourth order moment of the exact and approximate solutions differ at time t > 
of a quantity equal to B e t. Keeping the same notations of Subsection 14. 1[ we will use this fact to 
compute the ^3 distance between the scaled approximate solution h £ and the rescaled exact one h. 
Indeed, it is well know that if k € N, the k th moment of a distribution correspond to the k th derivative 
of its Fourier transform. Then g(-,t) and g £ (-,t) are at least four times differentiable and according 
to Taylor's theorem (|4.15|) - (|4.16|) imply 

Thus, one has for R > 



d 3 (h £ (t),h(t)) = sup |^|3 



(4.17) < A + ^ ^ R- 

y ' ~ R 3 24(1 + t) 2 

Optimizing (|4.17|) over R and gathering the result with inequality (|4.7|) (with s = 3) gives 
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Theorem 4.1. Let < go G L\ satisfy (|4.14|) . If g £ is solution to the associated Cauchy problem 
for Rosenau approximation \2. b]) with a kernel M £ verifying (|2.3|) . then the scaled distribution 

h £ (t, v ) ■.= FCtrvfV^rv) 

for V(t) := (1 + 1)~ 1 / 2 is such that, for all t > 

d 3 (M*W)< (^-^^3(20,^ ' , 

where uj a is the Gaussian distribution (|4.6p and 

B £ :=\l M £ (v)vUv. 

Remark 8 . The result of Theorem 14.11 indicates in a clear way that the rate of convergence towards 
the fundamental solution improves as soon as the initial value has a a higher number of moments 
equal to those of the Gaussian. In the present case, with four moments the difference in Fourier norm 
between solutions converges to zero at rate (1 + t)~ 3 / 4 . This suggests that the rate of convergence 
approaches the optimal rate when the number of finite moments of the initial value approaches 
infinity. This aspect of the problem is deeply connected to the analogous one studied in the central 
limit theorem, which is referred to as Berry-Esseen estimates. We remark that Berry-Esseen-like 
estimates for Fourier metric based metrics have been studied in |llj . 

Let us now specify the kernel M £ , to check if there are essential differences among the rate of 
convergence in e. 

The choice of the central differences kernel 

M £ {v) = ^[5 (v + ea) + S (v - ea)} 

gives B £ = 2e 2 <7 4 . Then 

1 / r \ 3/2 

d 3 (h £ (t),u a ) < (1 + t) 3 /2 ^(go,^) + Ci£ 3 / 2 lj^r t ) . 

where C\ is a positive constant depending on a. Compared to (|4.12|) . the rate in e is improved by a 
power 1/2. 

If we consider now the Rosenau kernel 

K(0 1 



i + M 2 e 2 ' 

B £ = Ke s for a nonnegative constant K. Then 



3/2 



;i + t)3/2^u,^, . \l + tj ' 

where C 2 is a positive constant depending on a. 

Remark 9. We can see here that the Rosenau kernel gives a better convergence rate in e towards 
the exact solution, compared to the central differences kernel. This improvement also depends of the 
number of finite moments. This shows that the number of finite moments plays an essential role also 
in connection with the order of convergence in e. Indeed in formulas (|4.12|) and (|4.13|) . which refer to 
the situation in which only moments up to order two have been considered, this difference between 
the approximations due to the two kernels was not evident. 
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4.3. Strong Convergence of a "Regularized" Approximate Solution. Let us assume that the 
kernel M £ belongs to The analysis of Section 0] shows that, in a suitable scaling which allows 

to maintain the energy of the solution bounded, there is convergence in Fourier distance towards a 
Gaussian function. This result, however, cannot be directly used to conclude that the solution to 
the Rosenau kinetic equation converges towards the fundamental solution to the heat equation, as 
time goes to infinity. In fact, since the d s -metric is not scaling invariant, and the decay in time of 
the distance found in Theorem 14. II is of order t -3 / 4 , by reverting to the original variables, the decay 
in time disappears. Nevertheless, these weak convergence estimates can be fruitfully employed to 
prove that the solution to Rosenau equation can be split into two parts, one of which is exponentially 
decaying to zero with respect to both time t and the parameter e, while the other converges strongly 
in time towards the fundamental solution to the heat equation at a lower rate. This result clarifies 
the nature of the Rosenau-type approximation to the heat equation. This approximation produces 
a solution which essentially consists of two parts, which differs for their regularity. The part with 
lower regularity (essentially the first term in (|3,4|) with the same regularity of the initial datum) is 
rapidly decaying to zero. The regular part is shown to behave like the classical solution to the heat 
equation, and decays towards the self-similar Gaussian solution for large times. The result will be 
achieved by making a heavy use of the representations given in Section O In particular, we will write 
the fundamental solution (|3.9|) by splitting it in two parts 



Xt/e 2 p M e (0\t/e 2 



1 - M, 



«>)) + ( 



l-M £ (0)e 



-Xt/e 2 



(4.18) 



£,reg 



-Xt/e 2 



As we shall see, this splitting separates in a natural way the singular part of the kernel from its 
regularized part P £ ,reg- 

Let g (respectively g £ ,reg) a solution to the heat equation (|1.4|) (resp. a solution to the Rosenau 
equation (|1.5jl obtained by convolution with the regularized kernel). In other words 

g(v,t) = g (v) * n<r (•,£)> 

(4.19) g £jreg (v,t) = g (v) * P E>reg (;t). 

We remark that, resorting to formula (|3.4j) . one recovers that g £treg (v,t) coincides with the Wild 
sum expansion in which the zero-order term is substituted by the regular term M e * go(v). Hence, 
9e,reg(v,t) is nothing but a regularized version of the Wild sum expansion (|3.4j) . Since the initial 
datum go(v) is of unit mass, 



\\g(t) ~ 9e,reg 

On the other hand, if A = a 2 

P £ , regit t)-h a (Z : t) 



(U a (v - w,t) 



e,reg 



(v - w,t)) g (w)dw 



dv 



< \\n a (t)-p £! 



< 



'■eg 



e ~o 2 t/e 2 ( e Ml(S)* 2 t/e 2 



e ~a 2 t/e 2 e M s {i)a 2 t/e 2 



(i-M £ (6V: ■ 



-a 2 \i\ 2 t 



+ 



1 - M £ (0 



-a 2 tje 2 



and this implies 

(4.20) d 2 (P £ , reg (t),Q a (t)) <d 2 (G £ (t),Q (7 (t)) + d 2 (M £ ,5 )e-° 2t / £2 . 

By means of the scaling property of the Fourier metric d s , formula (|4.2U|) ensures that, after the scaling 
f -> t^/VTTt, the convergence results (in scaled variables) of Section H] guarantee the convergence 
of the regularized kernel towards the Gaussian fundamental solution. 

Now, thanks to the regularity of P Streg (t), we can consider also convergence of the (regularized) 
approximate fundamental solution P ereg towards the heat kernel CT in stronger norms, typically L , 
which are invariant with respect to dilatation. Actually, it is enough to deal with I? norms. Indeed, 
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it is well known [HE] that if / G L 2 (R) is of finite kinetic energy, there exists an explicit, nonnegative 
constant C\ such that 



(4.21) 



L 1 



4/5 
L 2 



\v\ 2 \f{v)\dv 



1/5 



Note that the second moment of the quantities involved is uniformly bounded only in the scaled 
variables used in Section |U In original variables, in fact, the second moment of the solution, both of 
the heat equation and its Rosenau approximation, is increasing linearly with respect to time. Then, 
it is possible to interpolate the I? norm by some appropriate Fourier and homogeneous Sobolev 
norms. More precisely, it was shown in [3] that for all r G (0, 1) there exists an explicit, nonnegative 
constant C 2 = C 2 (r) such that 



(4.22) 



\\f-9\\i? <C 2 d 2 (f,g) 



2(1-0 



|2r 



\\f -9\\hn, 

where N = 3(1 — r)/r. Then, according to inequalities (|4.21|) - (|4.22|) - (|4.2U|) and to the results of 
Subsection 14. 2\ we will obtain the strong convergence of P £<reg (t, •) towards £l a (t, •) when t — ¥ 00 if 
this function belongs to the homogeneous Sobolev space H (R) for arbitrarily small N, and at the 
same time its Sobolev norm increases in time at a rate such that the right-hand side of inequality 
(|4.22p decays to zero as time goes to infinity. 

We will deal with the Rosenau kernel with A = a 2 

Kit) 1 



which belongs to H f 



1 + {eafe 

for all s < 1. Hence, in scaled variables, for all s < 1 



-a 2 t/e 2 



M, 



VT+t 



-u 2 t/e 2 



1+s 



sa 
2 



According to (|4.18|) . it remains to check the smoothness of 

gl^t) :=e ^)- 2 ^ 2 _l, V t >0. 
Let us set s < 1. We have thanks to de VHospital rule 



lim e +1 <3e{t,t) 



lim 



a 2 t 



e +2 M £ {£)e 



00 (s + l)e 2 

2aU e +3 



M e {i)a 2 t/e 2 



| ? |Too 8 + 1 (1 + (safe) 2 



a 2 t 



cxp 



e 2 (1 + (ea)H 2 ) 



0. 



Then we have proved that £ s {7 e (£;,i) = o(^ _1 ) which means according to Riemann criterion that 
G £ (-,t) G H S (R) for any t > 0. Detailed computations, reported in Appendix, allow to conclude that, 
for s < 1, 



(4.23) 



-a 2 t/e 2 



VT+t 



<C(e,s)(l + ty 



11 s 



where the explicitly computable constant C(e,s) depends on e, s and M e , and 5 > can be taken 
as small as we want. Thanks to inequalities (|4.20p - ([4.2ip - ([¥.22p - (|4.23p and Theorem 14.11 we obtain 
the following 

Theorem 4.2. Under the assumptions of Theorem \4-l[ if g is a solution to the heat equation (\1A\) 
and g e ,reg is given by (|4.19|) for the Rosenau kernel (|1.3p . then one has 



lim \\g(t) - g £ 

t— »oo 



reg 



0. 
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Remark 10 . Theorem 14.21 has various consequences. Coupling the decay of the d<i metric obtained 
in (|4.12|) with inequalities (|4.21|) and (|4.22|) . the regularity of P £jT eg implies a rate of decay of the L 1 - 
norm smaller or equal to (1 + t) 1 / 15 . On the other hand, the singular part of the fundamental solution 
t) converges to zero exponentially both in time and with respect to e. Hence, after a transitory 
time, there is prevalence of the regular part of the fundamental solution, which decays to zero in 
L l at a suboptimal rate. Hence, the main difference between the large-time behavior of the heat 
equation and its Rosenau approximation consists in a slower rate of convergence of the latter towards 
the Gaussian fundamental solution. The method of proof also indicates that the rate of convergence 
in L 1 is strictly linked to the regularity of the distribution of the background. For this reason, the 
argument leading to Theorem 14.21 can not be applied to the central difference approximation. Maybe 
a different approach will help to clarify if, in some weaker norm, a result similar to that of Theorem 
14.21 holds also for the central difference type approximation. 



5. Conclusions 

We studied in this article the validity of the approximation to the linear diffusion equation proposed 
by Rosenau [IB] as a regularized version of the Chapman-Enskog expansion of hydrodynamics. This 
approximation essentially is realized by substituting the heat equation with a linear kinetic equation 
of Boltzmann type, describing collisions of particles with a fixed background. This remark allows to 
consider the Rosenau approximation as a particular realization of a model Boltzmann equation, in 
which the background distribution is a general probability density with bounded variance. In addition 
to the Rosenau distribution, we considered in this article also a point masses background, which 
furnishes the central difference scheme to solve the linear diffusion equation. The main differences 
between the action of the two different backgrounds have been studied in some details. In particular, 
our analysis put into evidence that the approximation with a regular kernel is a good approximation, 
in that it possesses most of the typical properties of the heat equation, including the same large-time 
behavior, apart from a slower rate of decay towards the fundamental solution and a finite speed 
of propagation at high Fourier frequencies. We were not able to prove an analogous property for 
the non regular approximation, for which the rate of decay towards the equilibrium distribution has 
been only proven in scaled variables. The results of the present paper shall lead to new numerical 
approximations of the whole Rosenau approximation For instance, Wild sums have been 

successfully used in connection with the approximation of the bilinear Boltzmann equation (see e.g. 
|14j ) to develop a new family of Monte Carlo methods, which possess high order accuracy in time, 
and are asymptotic preserving (with respect to e). Following the same strategy, our next goal is to 
develop an efficient Monte Carlo method to compute the Rosenau approximation using our analysis 
on the linear Wild sums. 



Appendix A. Growth in Time of the Sobolev Norm of the Regularized Kernel 
For any given s < 1, we investigate the growth in time of the quantity 



y/T+t 



B s {t) = e~° l l £ 

To avoid inessential heavy notations, in what follows we fix e = a = 1. In this case 

t 



B s {t)=e 



-t 




+t)- i e 

By a simple change of variable into the integral, we can rewrite B s (t) in the equivalent form 

n2 ^ 1/2 



(A.i) B.(t) = {i + ty +l ' 2 i s (t) l i 2 = (i + ty +1 / 2 1 I \tf 



exp —t 



exp(— t) 
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The function I s (t) defined in (jA.ip satisfies the differential equation 
(A.2) d 
where 



I s (t) = -2I s (t) + 2 f A s (tt)d£, 
at Jr 



\e s ( , e 

^ exp [ -t 



exp —t 



exp(— t) 



For a given positive constant 5 < 2p + 1, let us set 



a = 1 



2p+ 1 



Then we find 
(A.3) 



If now |£| > (I + t) a we have 



e > (i+tr 2a 



i+e - 1 + (i + 1)- 2 " i + (i+t) 



2a 



>^(l + i)" 2a , 



and 



exp — i 



i + e 2 



exp(— t) 



exp - exp(-i) t 



< 



1 + ¥ ~ 1 + £ 



2 4 



The previous inequalities imply 

" + 00 



Id 2 * 



/ A s (e,t)(ie<texp(-t(I + t)- 2 72) / 7T r , 

(A.4) < c s texp(-t(I + t)- 2a /2) , 

where we denoted by c s the bounded constant 

|£|2a 

rTe 2 ) 2 ^' 

Putting together inequalities (|A.3|) and (|A.4j) we conclude with the upper bound 

(A- 5) := / A s (£,tK < o X i 

JgeR 2s + 1 (1 + t) zs+L 

Hence, considering that the second term in (|A.5p decay faster at infinity than the first one, for any 
given constant a > there exists a bounded constant C a such that if t > a 

(A.6) A s (t) < C a t-( 2s+1 - 5 \ 

Substituting inequality (|A.6[) into (|A,2j) . we obtain, for a given a > 0, and t > a 

I s (t) < /.(a) e~ 2t + e~ 2t [ A s (r) e 2r dr 



Mt) := J A s (£ ; t) d£ < 2s + l{ l + + 'L + c s t exp (-* (1 + r 2 72) • 



< J s(a ) e -2t + Ca e -2t /"* r -(2s + l-5) e 2r ^ 
J a 

f r^ 2s+1 - s h 2T dr. 

J a 



On the other hand, integration by parts shows that 

7a S 2 + 2a 

Choosing a in such a way that 

2s + 1 - S 1 

2a ~ 2' 
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we conclude that, for t > a 

e -2t [* T -{2s+l-5) e 2r dT < t -{2 a +l-S) i 
J a 

which implies 

(A.7) I s (t) < I s {a)e- 2t + C a t^ 2s+1 - & \ 

Using (|A.7j) into (|A.lj) gives 

(A.8) B s (t) = (1 + t) s+1 / 2 I s {t) 1 ' 2 < + t) s+1/ V + (1 + tf 2 VC~a- 

Since 5 is arbitrarily small, inequality (|4.23p follows. 
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